Quantum transport in chains with noisy off-diagonal couplings 
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We present a model for conductivity and energy diffusion in a linear chain described by a quadratic 
Hamiltonian with Gaussian noise. We show that when the correlation matrix is diagonal, the noise- 
averaged Liouville-von Neumann equation governing the time-evolution of the system reduces to the 
Lindblad equation with Hermitian Lindblad operators. We show that the noise-averaged density 
matrix for the system expectation values of the energy density and the number density satisfy 
discrete versions of the heat and diffusion equations. Transport coefficients are given in terms of 
model Hamiltonian parameters. We discuss conditions on the Hamiltonian under which the noise- 
averaged expectation value of the total energy remains constant. For chains placed between two 
heat reservoirs, the gradient of the energy density along the chain is linear. 



I. INTRODUCTION 

The problem of energy and charge transport in low- 
dimensional systems has attracted considerable current 
attention. 0,0, 0, S IE @| ■ This recent interest is pri- 
marily driven by new advances in nanotechnology and 
the resulting need to explain new and unusual transport 
properties of nanomaterials 0,0- I n many respects, such 
materials behave as effectively one or two dimensional 
systems. Furthermore, from a theoretical standpoint, 
this tnterest is driven by the desire to understand how 
dimensionality affects transport properties and whether 
normal transport is possible in low dimensions, and, if 
so, what are the conditions for such behavior. 

Studies of heat transport in one dimensional models 
have long and rich history. Review articles on heat 
transport in one dimension include Refs. 0] and [Tfj| . 
The earliest studies addressing this problem for anhar- 
monic chains include the seminal numerical simulations 
by Fermi, Pasta and Ulam E2 an d Peierls' discussion of 
the constraints on phonon interactions in one dimension 
fl2|. More recent theoretical studies of heat transport 
in one dimension include Refs. 

[3 13 13 El Gl Ex- 
tensive numerical simulations have also been reported in 
recent years [ToL UtI ITS IT^. |20| . It is now well established 
that thermal transport in momentum conserving one di- 
mensional models is anomalous, i.e. the Fourier heat law 
is not observed in this systems. The physical reason for 
such behavior is in the slow decay of normal modes with 
small wave numbers. In one dimension, the fast propa- 
gation of heat carried by these modes makes normal heat 
transport impossible. 

Due to the difficulties in the theoretical treatment of 
one-dimensional transport, exactly solvable models are 
of interest. One of the few one dimensional models that 
shows normal heat conductivity is a stochastic and non- 
Hamiltonian model introduced by Bolsterli et al. [2lLl2^| . 
In this paper we consider a quantum mechanical Hamil- 
tonian model in which the Hamiltonian matrix elements 
are taken as time dependent stochastic variables. Wc 
show that upon averaging over the noise, one obtains a 



Lindblad-type equation for the evolution of the quantum 
density matrix. We then go on to give the conditions 
upon the Hamiltonian which produce normal (i.e. diffu- 
sive) energy and particle transport. Finally, we demon- 
strate the validity of our model using numerical simula- 
tions. 



II. GAUSSIAN WHITE NOISE HAMILTONIANS 

The Hamiltonian of the one-dimensional model that 
we will consider in this paper belongs to a general class 
of Hamiltonians with Gaussian white noise 



H = Hn 



5>(i)^, 



(1) 



where Hq and Vi are arbitrary Hermitian operators and 
gaussian stochastic coefficients Xi(t) have average mean 
Xi(t) = and second moments given by 



\ i (t)\ j (t')=g ij 6(t-t'). 



(2) 



Hamiltonians such as this can model a wide variety 
of physical situations where the motion or transport is 
driven by an external field. One such example is the case 
of Forster resonant excitation transfer in biomolecules 
where the migration and diffusion of an initial electronic 
excitation within the system is dependent upon the lo- 
cal environmental and conformational fluctuations that 
modulate the off-diagonal couplings. 

From a theoretical standpoint, such Hamiltonians are 
useful since one can explicitly average over the noise 
when calculating the time evolution of the density ma- 
trix or expectation values of various observables. Noise- 
averaged time evolution for various specific forms of the 
Hamiltonian in Eq. (fj^) have been considered previously 
[3 13 HI El Ea- We give here a brief discussion of 
noise-averaged time evolution for the general case and 
then go on to emphasize points that are relevant for the 
studies of heat transport. 
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In general, the quantum density matrix satisfies the 
Liouville-von Neumann equation, 

l ^ = (C + C v (t))p. (3) 

Action of the superoperators Co and Cv on p is 
given, respectively, by Cop = j?[Ho,p] and Cy{t)p = 

I 



i J2i \ (t) [Viip]- The density matrix at time t is given 
in terms of the density matrix p(0) at time t = by 

p(t)=U(t)p(0). (4) 

Here the time-evolution superoperator U(t) is given by 
the following infinite series: 



U{t) 



-iCat 



i dTe- iCo ^ T, >C v (T)e- lCoT - dr dr 1 e-^-^ C v (T)e- iCo{ - T - T '^ C v {T')e- iCoT ' + 



(5) 



Noise-averaged expectation values of an operator O are 
computed using 

T0W)=^(0{t)U{t) P {Q)) . (6) 

It is assumed here that operator O can have an explicit 
time dependence. When performing averages as in Eq. 
© , we need to distinguish between two types of opera- 
tors: those with and those without stochastic coefficients 
\i{t). In the latter case, averaging over noise as in Eq. 
El reduces to an averaging of the evolution superoperator 
U(t) and such expectation values can be calculated with 
the noise averaged density matrix. 

Noise averaging of hi (t) can be performed by taking av- 
erages for each term in the series and then resumming the 
series. This involves averaging products of the stochastic 
coefficients. Since \i(t) is sampled from a Gaussian de- 
viate, all terms involving an odd number of coefficients 
necessarily vanish. Furthermore, any term with an even 
number of terms can be written as a sum of all possible 
products of second moments. However, due to the order 
of integrations over time in Eq. JSJ) and the fact that sec- 
ond moments in Eqs. @ involve delta functions in time, 
only one product from the sum contributes to the aver- 
age after all the time integrations are performed. For ex- 
ample, the fourth order average Ai(i)Aj(r)Afc(r')A;(r"), 
decomposes into the following second-order terms: 

A; (t)\j(T)X k (T')Xi(T") = g ijgM 5(t - t)5{t' - r") 
+g lk 9jiS(t - t')S(t - t") + g a g jk 5(t - t")6(t - r') 

(7) 

Since the region of integration for the fourth order term 
in the series © is t > r > r' > t" we can see that only 
the gijgki5{t — t)6(t' — t") term will contribute to the 
series. Similar analysis can be applied to all even order 
terms. 

Averaging over noise and resumming the series given 
by Eq. © produces 

Wj)=e~ lCot - M \ (8) 



where the action of the superoperator A4 on the density 
matrix p is given by 

Hence, our noise-averaged propagation is unitary and 
norm-conserving. It follows from Eq. JSJ that the noise- 
averaged density matrix satisfies the following equation 

= (Co - iM)p. (10) 

This allows one to see an interesting connection between 
noise-averaged time evolution of the density matrix for 
the noisy system and the time evolution of the reduced 
density matrix for an open quantum system. 

Consider the case when correlation matrix g^ in Eq. 
© is diagonal, i.e. <7y = guSij. In this case, Eq. ltTU|) 
can be rewritten as 

i 

This has the form of a Lindblad equation with Hermitian 
Lindblad operators |2fi|. Thus, the time evolution of the 
reduced density matrix under the Lindblad equation with 
Hermitian operators is equivalent to noise-averaged uni- 
tary evolution with Hamiltonian Q in the special case 
of a diagonal correlation matrix g^ . Connection between 
Lindblad-type evolution with Hermitian Lindblad opera- 
tors and an averaged random unitary evolution was dis- 
cussed in Ref. [2!j. It is interesting that the unitary 
operator proposed in Ref. [2fll | is different from the uni- 
tary operator corresponding to the Hamiltonian in Eq. 
Q). This implies that averaging different unitary evolu- 
tions over the same noise spectrum can ultimately lead 
to the same Lindblad-type evolution. 

If the operator whose noise-averaged expectation value 
we want to find depends on the coefficients Xi(t) then 
these coefficients must be taken into account when per- 
forming averages as in Eq. [S] Let us consider the linear 
case Xi(t)W where W is time independent and Xiit) is a 
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gaussian random variable. An example of an operator of 
this type is the second term of Eq. ^ Using series (JSJ 
for superoperator U{t) and substituting \{t)W for O in 



Eq. (|SJ| we can again perform averages for each term and 
resum the series to obtain 



(12) 



where the last line is obtained by using cyclic permutations under the trace and the fact that e lCot Mt p(0) — p(t). 



Explicit averaging over noise was performed in Eqs. 
JSJ) and i|12|) . However, these expressions remain mostly 
formal. For some forms of operators H and Vi further 
exact treatment is possible. One such example is Hq and 
Vj, that are quadratic forms of creations and annihilation 
operators. Hamiltonians of this type were used in solid 
state physics to describe electronic transport under the 
influence of an external driving field such as lattice vibra- 
tions [2(| |3(j • These models gave a satisfactory picture 
of electron diffusion. However, such models generally do 
not conserve average energy of the system and therefore 
cannot show normal heat propagation. 

Let us use Eqs. (|SJ) and i|12[) to determine the precise 
conditions that the operators Hq and Vi must satisfy in 
order for the noise-averaged expectation value of the total 
energy to remain constant. The noise- averaged energy 
expectation value is given by 



{H(t)) = {H (t)) + j2mm)) 



(13) 



Since Hq does not depend upon the stochastic coeffi- 
cients, we obtain 



(H (t)) = Tr(H e 



-iCnt-Mt 



P(O))- 



(14) 



To find the second term on the right hand side of Eq. 
(|T3)l we use Eq. $T2$ where we need to put Vi for W and 
sum over all i. We obtain 



HtW(t)) = -~Tr p(t) Y,9ii [Vi, Vi\ (15) 



Using the fact that the trace is invariant under cyclic 
permutations, the symmetric nature of matrix g^j, and 
that CqHq — 0, we obtain 



d(H (t)) 
dt 



= -Tr lp(t)MHt 



(18) 



For arbitrary p(0) this derivative will be equal to zero if 
MH ° = ^2 E 9ij [Vi, [Vi,H ]] - 0. (19) 



This is the condition for (H(t)) to be constant in time. 
Examples of operators Hq and Vi for which Eq. I|19(l is 
true include the case when Hq and V, satisfy the condi- 
tion [Hq, Vi] = const for all i, in particular when Hq and 
V. commute. 



III. MODEL HAMILTONIAN AND 
OBSERVABLES 

We now apply the results of the previous section to 
study quantum heat conduction and diffusion in a one- 
dimensional chain of N sites described by the following 
Hamiltonian 



JV 



H = ^y^qtaj + %y^ j p, i (t)(a\a i+ i + a\ +1 a,i) 

i i—1 
N 

+ifi Vj (t) (a\ai + i - a\ +1 ai) (20) 



Since matrix gij is always symmetric, the last quantity 
vanishes f or arbit rary V^'s and, therefore, does not con- 
tribute to (H(t)). Thus we have 



(H(t)) - (H (t)) 



(16) 



For (H(t)) to be constant in time, its time derivative 
must vanish. Differentiating (ffo(i)) with respect to time 
we obtain 



d(H (t)) 



dt 



Tr(H (-iC -M)p(t)). 



(17) 



Here a\ and denote, respectively, creation and annihi- 
lation operators for quasi-particles at chain site i. These 
quasi-particles can be bosons, fermions, or spin excita- 
tions. The stochastic hopping coefficients p,i(t) and Vi{t) 
are realizations of the Gaussian noises with zero mean. 
We choose the following simple expressions for the second 
moments, 



H(t)p k (t>) = aiS jk 5(t-t'), 



p,j(t)i>k(t') 



a 2 2 5 jk 8{t - 1') 




(21) 
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We will consider two types of boundary conditions: cyclic 
and open ends. For the former we have ajy+i = a±, 



t 

Af+l 



a[, and for the latter ajy+i 



a 



N+l 



0. 



The model Hamiltonian in Eq. I|2UI) belongs to the 
class of noisy Hamiltonians we just described. Indeed, 
we can rewrite Eq. l2*Ulas 



27V 



(22) 



where 



N 



a]a,i 



and 



Xj(t)Vj = hm(t)(4a 



i+1 



l a i) 



for even j's and 

x j(t)Vj = ihvi(t)(a\ai + i - a\ +1 ai) 

for odd j's. Thus, discussion of Sec. [H]can be applied to 
the studies of time evolution driven by this Hamiltonian. 
Note that Hamiltonian l|20|) satisfies the average energy 
constraint in Eq. Ijl9|l because its time independent part 
Hq commutes with all Vi's. Therefore noise averged ex- 
pectation value of the total energy remains constant in 
time and we can discuss heat conduction in this model. 

Heat conduction and particle diffusion can be studied 
by considering time evolution of expectation values of 
energy density and number density averaged over noise. 
For both types of boundary conditions, operators for the 
number density rii are 



(23) 



The energy density ti at site i for cyclic boundary con- 
ditions is given 



hujajai + ^%Hi(t)Ai + -7i/i,_i(t)Aj_i 



(24) 



where Ai = {a\a i+ i+a\ +l a,i) and B t = i{a\a l+ i-a\ +1 ai) 
In the case of open-end boundaries, operators for the en- 
ergy density at site i are the same as for cyclic boundaries 
except for the first and last sites where we use 

€\ = huia\ai + —h/ii(t)Ai + —hi'i(t)Bi, 

e N = huja^ N a N + + ^hv N -i(t)B N -i 

(25) 



noise independent component (which is proportional to 
number density) and a noisy component linear in coef- 
ficients fii(t) or Vi(t). Thus, we will need to take both 
typ es of av erages over noise discussed in Sec. |^ We have 
for (ni(t)), 



(m(t)) = Tr(nie 



-iC t-Mt 



(26) 



Following procedure in Sec. ^and using the explicit form 
of the second moments (Eq. I21|) we obtain, 

Mr = f^[4 [A u p]] + Pi'PU ■ ( 2? ) 



Here the summations extend from 1 to N for cyclic 
boundaries and from 1 to N — 1 for open-end boundaries. 
In the case of (ei(t)), we can evaluate the noise-dependent 
part using Eq. Substituting Ai or Bi for W and using 
Eq. |^we obtain 



{lM(t)Ai(t)) = {"i{t)Bi{t)) = 0. (28) 

Thus, for the special form of second moments in Eq. 121|l , 
time dependent part of does not contribute to (ei(t)) 
and we have 



i(t)> = hu(m(t)). 



(29) 



In order to arrive at equations of motion for {ti{t)) and 
(rii(t)), we differentiate both sides of Eq. (|26|l with re- 
spect to time and obtain 



Of 



= -Tr(ni(iCo + M)e-* Cot - Mt p{0)) 

= -Tr(pJt)Mni). (30) 



To obtain the last expression we used explicit form of su- 
peroperators Co and M, cyclic permutations under trace 
and the fact that Corii — 0. Using the explicit form of 
superoperator Ai we obtain the following equations for 
M*)) and (e;(i)) 



d(m(t)) 



dt 



a z ^{n % {t)) 



d(ei(t)) 



dt 



ofAiieiit)) 



(31) 
(32) 



here a 2 = a\ + a\ and is a discrete Laplacian whose 
action on some function fi defined over discrete set of 
points i is given by 



&ifi = -2fi + fi+l+fi-i. 



(33) 



In the case of cyclic boundary conditions Eqs. I3H and 1321 
are valid for all i's, whereas for open end boundaries they 
are valid for all i's except i = 1 and i — N for which we 
have 



Note that the number density operator is independent 
of the noise. The energy density, however, has both a 
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d(n N (t)) 
dt 



ar((ntf_i(i)) - (njv(i)», 



dt 
0t 



= ^((e a (t)>-(ei(t)», 



= ^(( ew _!(t)> - Mt)» 



(34) 



A, is the discrete analogue of the continuous operator 
l 2 d 2 /dx 2 where I is interatomic distance. We can see that 
for cyclic boundary conditions, (ftj(t)) and (ej(i)) satisfy 
discrete versions of the diffusion and heat equations, re- 
spectively. For open end boundary conditions these equa- 
tions are satisfied everywhere except the boundaries. In 
both cases the stationary solutions correspond to con- 
stant (rij(t)) and (ej(t)). Coefficients of diffusion and 
energy density diffusion are given by 



K d = K h = a 2 l 



2,2 



(35) 



In the case of diffusion coefficient Kd, Eq. 13511 agrees 
with the earlier results |U Time dependent dif- 

fusion coefficient obtained in Ref. [24j reduces to our 
result Eq. (|35|l for the special case of Hamiltonian l|20() 
because all time dependent terms vanish for Hamiltoni- 
ans with zero averaged hopping coefficient. We also note 
here that the closed equations for (rij(t)) and (ei(t)) are 
due to the simple form of Hamiltonian (|20() and the re- 
sulting structure of superoperator A4. For more general 
Hamiltonians, such as those considered in Refs. jUEl], 
closed equations exist only for the noise-averaged single 

particle density matrix, (aj<2j). However, such more gen- 
eral models do not conserve noise averaged energy of the 
system and therefore cannot be used for studies of heat 
conductivity which is the primary focus of this paper. 
Including on-site noise terms in Hamiltonian i|20|) does 
conserve the noise averaged energy. However, inclusion 
of these terms does not affect the transport equations 
given by Eqs. JSH EH»- 



IV. COUPLING TO THE RESERVOIRS 



sites to Hamiltonian l|20|) and, secondly, by introducing 
a non-Hamiltonian damping for the creation and anni- 
hilation operators at the terminal sites with decay rates 
of, respectively, 71 and 7 at. The noise sources ?7i(i) and 
r?Ar(i) are assumed to be not correlated to other noises in 
the Hamiltonian (|20|l or to each other and hence can be 
treated on the same basis as the hopping integral noises. 
We assume that 



vi(t)vi{t') = am - f), m (t) m (t') = e N s(t - 0.(36) 

Note that the stochastic coupling that we consider cor- 
responds to reservoirs that exchange both energy and 
particles with the lattice. In the case of bosons, modified 
equations for energy density of the first and last sites are 
given by 



d(ei) 
dt 

dt 



= -(2 7l + a 2 )( ei )+a 2 (e 2 )+^ 



(37) 



= -{2 lN + a 2 )(e N ) +a 2 (e N - 1 ) + ^ (38) 



Equations of motion for energy densities at other sites 
remain unchanged and are given in Eq. I|32|l . The sta- 
tionary solutions of equations correspond to a linear dis- 
tribution of (ej). The energy gradient 9 



M - (ei) 



N -1 ' 

is given in terms of reservoir coupling parameters by 



(39) 



2(N - 1)71 7at + a 2 (71 + 7jv) 
while for the stationary state 



7-r w£ + o?t 



271 



N 



a 2 9 



271 



(40) 



(41) 



In the case of fermionic or spin operators, equations of 
motion for (ei) and (ejv) are 



We would like to verify if a linear energy density dis- 
tribution will be achieved if the lattice is placed between 
two stochastic Langevin-type reservoirs. To model cou- 
pling of the lattice to such reservoirs we will use our 
original open-end lattice model and modify the equa- 
tions of motions for the operators for the first and last 
sites. These modifications are achieved by, firstly, includ- 
ing additional time-dependent noise terms of the form 
hrji(t)(a\ +ai) and Kr]N(t)(a\f +ajy) for the first and last 



d(ei) 
dt 
dj^) 
dt 



-(2 7l + a 2 + 2lu£)(e 1 ) + a 2 (e 2 ) + ujg 



-{2 lN + a 2 + 2w£%)(e N ) + a 2 {e N ^) + u){%. 

(42) 



Again, the s tatio nary state corresponds to the linear dis- 
tribution of (ei) with the energy density gradient given 
by 



w(7i& 



2{N - l)( 7l - ugX-lN + w&) + « 2 (- 



-7i - In + u£ 2 + wff ) 



(43) 
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and 



a 2 9 



2(7i + ) ' 



a 2 



2( 7JV + ^) ' 



(44) 



Since in this case we are dealing with non-equlilibrium 
situation there is no unique way of introducing temper- 
ature in this case. However, we can say that the local 
temperature is proportional to the average energy den- 
sity at a given site. Hence, we recover the macroscopic 
heat equation starting from a microscopic quantum de- 
scription. 

V. DISCUSSION 

In this paper, we developed a quantum mechanical 
model for energy and particle transport in a linear chain 
with stochastic hopping coefficients. We present a set 
of specific conditions in which the noise-averaged energy 
expectation value will remain constant. Under such con- 
ditions we arrive at a quantum mechanical law for heat 
transport. Furthermore, we consider the case in which 
the chain can exchange both particles and energy with 
reservoirs on either end. Both the energy density gra- 
dient and particle density gradient can be obtained in 



closed form. We consider only gaussian noise for our 
stochastic coefficients. More detailed analysis will be re- 
quired to develop a complete understanding of how well 
the averages represent the behavior of individual mem- 
bers of the ensemble. To reconfirm our theoretical re- 
sults we performed a series of numerical simulations for 
our model l)20ll with a particle placed at one of the sites 
at t = 0. As expected normal diffusion was observed 
for noise averaged time evolution for wide range of noise 
strengths. The simplicity of this model allows us to de- 
velop a more general treatment of quantum transport in 
a wide range of physical systems ranging from resonant 
electronic energy transport in DNA chains and organic 
glasses to Joule heating in molecular wires. 
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